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ABSTRACT 

We present the X-shooter Lens Survey (XLENS) data. The main goal of XLENS 
is to disentangle the stellar and dark matter content of massive early-type galaxies 
(ETGs), through combined strong gravitational lensing, dynamics and spectroscopic 
stellar population studies. The sample consists of 11 lens galaxies covering the red- 
shift range from 0.1 to 0.45 and having stellar velocity dispersions between 250 and 
380 km s _1 . All galaxies have multi-band, high-quality HST imaging. We have obtained 
long-slit spectra of the lens galaxies with X-shooter on the VLT. We are able to dis¬ 
entangle the dark and luminous mass components by combining lensing and extended 
kinematics data-sets, and we are also able to precisely constrain stellar mass-to-light 
ratios and infer the value of the low-mass cut-off of the IMF, by adding spectroscopic 
stellar population information. Our goal is to correlate these IMF parameters with 
ETG masses and investigate the relation between baryonic and non-baryonic matter 
during the mass assembly and structure formation processes. In this paper we pro¬ 
vide an overview of the survey, highlighting its scientific motivations, main goals and 
techniques. We present the current sample, briefly describing the data reduction and 
analysis process, and we present the first results on spatially resolved kinematics. 

Key words: dark matter - galaxies: ellipticals and lenticular, cD - gravitational 
lensing:strong - galaxies: kinematics and dynamics - galaxies: structure - galaxies: 
formation 


1 INTRODUCTION 

Massive early-type galaxies (ETGs, cr Js 200 km s“ ) contain 
more than half of the stellar mass in the Universe (e.g., Bell 
et al. 2003). Their internal structure and dynamics as well 
as the relationship between baryonic luminous matter and 
dark matter therefore provide key quantitative tests of the 
ACDM cosmological paradigm (Blumenthal et al. 1984) and 
are crucial to fully comprehend the processes involved in the 
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hierarchical formation of galaxies (e.g. White & Rees 1978; 
Davis et al. 1985; Frenk et al. 1985). 

In hierarchical structure formation models massive 
galaxies are believed to form through mergers of lower-mass 
galaxies (e.g., Blumenthal et al. 1984), however their stel¬ 
lar populations are old and evolved (see Renzini 2006 and 
Conroy 2013 for recent reviews). In addition, analysis of 
the evolution of the stellar mass and the luminosity func¬ 
tions has demonstrated that there is little evolution in the 
co-moving space density of massive passive galaxies since 
2 ~ 1 (Cimatti et al. 2004; McCarthy et al. 2004; Glazebrook 
et al. 2004; Daddi et al. 2005; Saracco et al. 2005; Bundy, 
Treu & Ellis 2007; Perez-Gonzalez et al. 2008; Marchesini 
et al. 2009; Ferreras et al. 2009; Carollo et al. 2013; Muzzin 
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et al. 2013; Lundgren et al. 2014). The prevailing idea is 
that ETGs formed the bulk of their stars early in the evolu¬ 
tion of the Universe (i.e., cluster ETGs z > 3, field ETGs at 
2 > 1.5 — 2, Thomas et al. 2005; Renzini 2006; Cimatti et al. 
2008; Whitaker et al. 2013) and then only at lower redshift 
(i.e., 2 < 1 ) may have merged together to build the massive 
ETGs that we see today in the nearby Universe (Johansson, 
Naab & Ostriker 2012; Nipoti et al. 2012; Buitrago et al. 
2013). 

To reconcile these apparently conflicting results and 
to paint a more robust physical picture for the formation 
and evolution of these massive systems, enormous effort has 
been expended in the last two decades to constrain the rel¬ 
ative contributions of baryonic, dark matter and black hole 
constituents of ETGs through stellar dynamical tracers, X- 
ray studies, and gravitational lensing (e.g. Fabbiano 1989; 
Mould et al. 1990; Saglia, Bertin & Stiavelli 1992; Franx, van 
Gorkom & de Zeeuw 1994; Carollo et al. 1995; Arnaboldi 
et al. 1996; Rix et al. 1997; Loewenstein & White 1999; 
Gerhard et al. 2001; Romanowsky et al. 2003; Treu & Koop- 
mans 2004; Treu et al. 2006; Cappellari et al. 2006; Thomas 
et al. 2007; Czoske et al. 2008; Auger et al. 2010a,b; Coccato 
et al. 2010; Barnabe et al. 2011; Das et al. 2011; Treu 2010; 
Courteau et al. 2013). A variety of recent observations sug¬ 
gest that baryonic luminous matter, which dominates astro- 
physical observables, and dark matter, which dominates the 
dynamics during galaxy formation, are strongly linked. As 
a consequence, progress in understanding ETGs and more 
generally the theory of galaxy formation and evolution re¬ 
quires the precise measurement of their stellar and dark 
matter density profiles, together with measurements of how 
these quantities depend upon properties such as the mass 
and the age of the stellar assembly. 

Lensing and dynamics observations show that in the in¬ 
ner region of massive ETGs (inside few effective radii, J? e ff), 
where baryonic and dark matter are both present, the total 
mass density profile is well described by a power-law den¬ 
sity profile, ptot oc r 1 , with an almost isothermal slope 
of 7 ' « 2 and ~ 10% intrinsic scatter (e.g. Gerhard et al. 
2001; Treu & Koopmans 2004; Koopmans et al. 2006, 2009; 
Coccato et al. 2009; Auger et al. 2010a; Barnabe et al. 2009, 
2011; Napolitano et al. 2011; Bolton et al. 2012; Sonnenfeld 
et al. 2012; Dutton & Treu 2014). However, in the same re¬ 
gion, the dark matter density profile seems to be shallower, 
consistent with a density slope 1.0 < 7 dm < 1-7, although 
this slope is less well constrained than the total mass den¬ 
sity slope (e.g. Treu & Koopmans 2004; Dye & Warren 2005; 
Grillo 2012; Sonnenfeld et al. 2012, 2014, 2015). 

In addition, recent observations as well as theoretical 
studies based on stellar population and dynamical models 
(e.g. Bullock et al. 2001; Padmanabhan et al. 2004) indicate 
that, for a universal IMF , the dark matter fraction (/dm) in 
the internal region increases monotonically with the mass of 
the galaxy (e.g. Zaritsky, Gonzalez & Zabludoff 2006; Auger 
et al. 2010 b), a trend that is more conspicuous in the case 
of slow-rotator ellipticals (Tortora et al. 2009). On the other 
hand, lensing and dynamics studies also suggest that the 
luminous stellar mass-to-light ratio (T*) scales with the lu¬ 
minous mass of the system (Davis et al. 1985; Bardeen et al. 
1986; Bell & de Jong 2001; Girardi et al. 2002; Napolitano 
et al. 2005; Grillo et al. 2009; Auger et al. 2010a), under the 


assumption of standard Navarro-Frenk-White dark matter 
halo profile (NFW, Navarro, Frenk & White 1996). 

Despite this progress, it remains difficult to separate 
the stellar and dark matter components, mostly due to a 
relatively poor understanding of the precise shape of the 
stellar IMF and its associated Y*. Uncertainties related to 
the latter can easily lead to a factor of two uncertainty on 
the inferred stellar mass. 

Although it is commonly assumed that the IMF is uni¬ 
versal and independent of cosmic time (e.g., Kroupa 2001; 
Chabrier 2003; Bastian, Covey & Meyer 2010), several au¬ 
thors have suggested that the IMF might evolve (Dave 2008; 
van Dokkum 2008) or depend on the stellar mass of the 
system (e.g. Worthey 1992; Trager et al. 2000; Treu et al. 
2010; Graves, Faber & Schiavon 2009; Graves & Faber 2010; 
Auger et al. 2010b; Napolitano, Romanowsky & Tortora 
2010; Cappellari et al. 2012; Martfn-Navarro et al. 2015; 
Posacki et al. 2015). Recent results based on single stel¬ 
lar population (SSP) modelling of galaxy spectra, indicate 
that the number of low-mass stars (M < 0.3 Mq) increases 
more rapidly with galaxy mass than the number of high- 
mass stars (van Dokkum & Conroy 2010, hereafter vDCIO, 
Spinicllo et al. 2011, 2012, 2014; Cappellari et al. 2012; Tor¬ 
tora, Romanowsky & Napolitano 2013; Ferreras et al. 2013; 
La Barbera et al. 2013). This could imply that part of the 
increase in Y* with galaxy mass can be due to a changing 
stellar IMF rather than only an increasing dark matter frac¬ 
tion, consistent with the earlier findings by Treu (2010) and 
Auger et al. (2010b), based on combing lensing, dynamics 
and stellar population information. The shape of the IMF 
and the internal fraction of dark matter mass inside massive 
ETGs can not be both universal at the same time, but it 
remains difficult to break the degeneracy between these two 
effects. 

The X-shooter Lens Survey (XLENS, Spiniello et al. 
2011 , 2012 ) aims to take these analyses one step further. 
With a combined analysis of strong gravitational lensing, 
dynamical and spectroscopic stellar population studies, the 
XLENS project will be able to disentangle the stellar and 
dark matter content of galaxies and, for the first time, di¬ 
rectly constrain the normalization, shape and cut-off mass 
(Afi ow ) of the low-mass end of the IMF and correlate these 
results with other galaxy properties such as galaxy mass, 
size, stellar density and/or stellar velocity dispersion (e.g., 
Spinicllo et al. 2011, 2012, 2014; Barnabe et al. 2013; 
Spinicllo et al. 2015). 

In particular, our methodology consists of combining 
spatially resolved kinematics with high-precision strong- 
gravitational-lensing measurements of the total mass and 
the mass density profile near the lens Einstein radius (i?Ein), 
to obtain a precise internal dark-matter mass fraction in¬ 
side one effective radius. By taking advantage of the wide 
wavelength coverage and throughput of the X-shooter spec¬ 
trograph, we are able to obtain high-resolution spectra from 
the UVB to the near-IR for detailed stellar population and 
spatially-resolved kinematics studies of the high-mass end of 
the early-type galaxies (ETGs) up to redshift 2 ~ 0.7. 

A previously under-appreciated result of our approach 
is that it allows us to tackle the long-standing problem of 
constraining the low-mass cut-off (Mi ow ) for the IMF. In 
previous studies, the value of Mi ow has always been treated 
as a relatively unconstrained parameter, despite being crit- 
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ic.al to determine the value of T*. Stars with masses be¬ 
low ~ O.15A/0 have very little effect on the spectral lines 
and on the line-index measurements in the optical and near- 
infrared for any assumed IMF slope (Conroy & van Dokkum 

2012) , but they provide a substantial contribution to the to¬ 
tal mass budget of the system for IMF slopes steeper than 
~ 2 (Worthey 1994). Therefore, to constrain the value of 
Mow it is necessary to support spectroscopic studies with a 
combined lensing and dynamics analysis (see Barnabe et al. 

2013) . 

In this paper we provide an overview of the survey, 
focusing on the data reduction and analysis process. We 
present scientific results on the IMF slope and cut-off mass 
of the full sample in Spiniello et al. (2015) and other two 
companion papers of the Survey, in preparation. Moreover, 
we already presented results on the IMF slope of a single 
very massive lens galaxy used as pilot program to test our 
method in the XLENS I, Spiniello et al. (2011). In Section 2 
we introduce the XLENS Survey, highlighting its main goals 
and characteristics and we present its selection criteria. In 
Section 3 we describe the observations and data reduction 
process. In Section 4 we present the spatially resolved kine¬ 
matics up to ~ 1 R e s obtained from VLT X-shooter spectra 
of 10 ETGs with <r* 250kms _1 . We summarize our find¬ 

ings and present our preliminary conclusions in Section 5. 


2 OVERVIEW OF THE XLENS SURVEY AND 

TARGET SELECTION 

The XLENS aims to study the stellar population of strong 
early-type lens galaxies, for which exquisite HST multi-band 
images and very detailed lens models are available. The main 
goal of the survey is to disentangle the dark matter and stel¬ 
lar mass distributions, as well as to constrain the slope and 
the cut-off mass of the low-mass end of the stellar IMF di¬ 
rectly from spectra, by combining the lensing and dynamical 
results with the fully independent inferences from spectro¬ 
scopic stellar population studies. 

Our current sample has been selected from the Sloan 
Lens ACS Survey (SLAGS, Bolton et al. 2006, 2008a,b; Treu 
et al. 2006, 2009; Koopmans et al. 2006; Gavazzi et al. 2007, 
2008; Auger et al. 2009, 2010b; Newton et al. 2011). SLACS 
is a survey of lenses spectroscopically selected from the Sloan 
Digital Sky Survey (SDSS, York et al. 2000) with follow-up 
imaging conducted with the Advanced Camera and Spec¬ 
trograph (ACS) aboard the Hubble Space Telescope (HST). 
SLACS was initiated in 2003 with the purpose of confirming 
and imaging galaxies that act as strong gravitational lenses 
of emission-line background sources. The survey has greatly 
expanded the number of strong gravitational lens galaxies 
at redshift 2 ^ 0.5 with complete redshift information, pro¬ 
viding an ideal sample of systems to be analysed with joint 
lensing and dynamics techniques. In particular, half of the 
lenses were selected from the SDSS luminous red galaxy 
(LRG) spectroscopic sample (Eisenstein et al. 2001). The 
LRG sample is defined by photometric selection cuts that 
very efficiently select massive ETGs in the redshift range 
0.15 < 2 ^ 0.5, as confirmed by SDSS spectroscopy. These 
galaxies are very homogeneous in their spectral, photomet¬ 
ric, and morphological properties. The remaining targets 
were selected with the same spectroscopic selection proce¬ 


dure from within the MAIN galaxy sample of SDSS (Strauss 
et al. 2002). The MAIN sample is more heterogeneous, and 
therefore a quiescent, absorption-dominated spectral crite¬ 
rion was imposed by requiring the lens candidates to have 
rest-frame equivalent widths in Ha of EWh 0 < 1-5 A. The 
method by which lens candidate were selected is extensivey 
described by Bolton et al. (2006). 

The current SLACS sample consists of 131 strong grav¬ 
itational lens candidates (Bolton et al. 2008a; Treu et al. 
2009; Auger et al. 2009; Shu et al. 2015). To date, this 
represents the largest uniform sample of strong gravita¬ 
tional lenses and it includes 98 “grade-A” and 33 “grade¬ 
s’’ and “grade-C” (likely lenses) confirmed strong galaxy- 
galaxy lens systems complete with lens and source redshifts, 
F814W lens-galaxy photometry, gravitational lens models, 
and measured stellar velocity dispersions (Auger et al. 2009). 
Approximately 80% of the “grade-A” systems have elliptical 
morphologies while ~ 10% show some spiral structure; the 
remaining lenses have lenticular morphologies. Such a large 
and high-quality lens sample provides a unique resource for 
the quantitative study of massive galaxies. 

SLACS is a lensing-selected sample and hence also 
mass-selected sample, largely peaking around 1 — 2M* galax¬ 
ies. Shu et al. (2015), published after the XLENS data was 
acquired, increased the final sample extending it to lenses 
below below M* 

From the final SLACS sample of ~ 100 confirmed lenses, 
the XLENS sample includes only “grade-A” systems with el¬ 
liptical morphology that he above the knee of the luminosity 
function (<t*,sdss > 250 kms -1 ) 1 . 

The velocity dispersion cut has been made because of 
the finding that there is a trend of increasing fraction of 
“non-luminous” mass beyond M* galaxies (Treu et al. 2010; 
Auger et al. 2009, 2010a; Barnabe et al. 2013). This mono- 
tonically increasing of the mass-to-light ratio as a function 
of the velocity dispersion could be due to an increasing num¬ 
ber of low-mass stars (i.e. a low-mass IMF steeper than the 
Milky Way, Kroupa 2001; Chabrier 2003) or to an increase 
in internal dark matter, or a combination thereof. 

We note that this fraction is very low (less than 10%) 
around M* (Auger et al. 2009, 2010a). Although below M* 
this fraction might go up again, as suggested by other dy¬ 
namical and gravitational lensing measurements, such low- 
mass lenses are rare and much harder to study (see Shu 
et al. 2015). The current XLENS sample does therefore not 
include any of these only recently discovered lower-mass lens 
galaxies. 

The current XLENS-SLACS sample consists of eleven 
luminous red galaxies (including one previously published 
system not selected from SLACS) with apparent V- 
magnitudes between 15.8 and 18.1 (from Auger et al. 2009), 
high-resolution and high signal-to-noise (S/N) HST imaging 
in B, V, I and H-bands, and with detailed lensing models 
constructed by the SLACS collaboration. 

In addition to the lens galaxies selected from SLACS, 
another system has been observed in X-shooter Guaranteed 
Time (GTO, P087.A-0620): SDSSJ1148+1930, also known 
as the “Cosmic Horseshoe” (Belokurov et al. 2007). This 


1 SLACS is a lensing-selected sample and hence also mass- 
selected sample, largely peaking around 1 — 2M* galaxies 
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Figure 1. Parameter space of the XLENS sample (XLENS from 
SLAC: black triangles, Pilot Program: black star) compared with 
other ETG samples (colored box). The SLACS systems are plot¬ 
ted as red crosses. Some overlap at the high-mass exists between 
all these studies but the XLENS galaxies extend the redshift range 
until a lookback time of 5 Gyr. 

lens at redshift 2 = 0.444 was analysed as a “pilot pro¬ 
gram” for the survey. Results on this extremely massive 
galaxy ((<T*)(5i-Reff) = 352 ± 10 ± 16kms _ ) were presented 
by Spiniello et al. (2011). In the following sections, we limit 
our discussion to the ten SLACS-selected XLENS systems 
(with 0.11 ^ a ^ 0.44). 

In Figure 1 the mass-redshift space covered by the fi¬ 
nal XLENS sample (black triangles and star) is compared to 
those from the SLACS parent sample (red crosses) as well as 
other sample of massive ETGs that have been used to study 
the low-mass end of the IMF (Cappellari et al. 2013; Con¬ 
roy & van Dokkum 2012; La Barbera et al. 2013; Ma et al. 
2014, colored boxes). The different samples overlap in the 
high-mass regime but XLENS pushes the redshift bound¬ 
aries beyond the local Universe. 

The current sample of lens ETGs with their photometric 
parameters and physical properties is presented in Table 1. 


3 OBSERVATIONS AND DATA REDUCTION 

X-shooter (Vernet et al. 2011) observations of the all systems 
were carried out in a total of 30 hours of Observation Time 
(OT), awarded in two periods (P086 and P089) 2 . The long- 
slit mode was used, splitting the beam over three arms: UVB 
(R=3300 with a 6 slit); VIS (R=5400 with a l”5 slit); and 
NIR (R=3300 with a l'/5 slit), covering the wavelength range 
from 3000 to 25000 A simultaneously. The 11" long slit was 
always centred on the lens galaxy with the position angle 
(PA) aligned with the major axis of the systems or with 
an angle that minimizes the contamination from the lensed 
background source and leaves enough sky region to facilitate 
accurate sky subtraction. A different number of Observation 
Blocks (OBs) have been obtained for each system to reach 

2 15 hours OT during P86 (P086.A-0312, PI: Koopmans) and 15 
hours OT during P89 (P089.A-0364, PI: Spiniello). 


a S/N high enough to perform stellar population analyses 
(S/N ~ 50 per A). Each OB had the same layout with three 
scientific exposures for a total exposure time on target for 
each arm of ~ 2500 s per OB 3 . During the observations, 
the seeing varied from ~ Ql!Q to ~ Standard calibration 
frames and standard stars for flux calibration were obtained 
after each OB. 

The data have been reduced using the ESO X-shooter 
pipeline v. 2.0.0 (Goldoni et al. 2006) and the Gasgano 
data file organiser developed by ESO. The pipeline reduction 
performs standard bias subtraction and flat-fielding of the 
raw spectra. Cosmic rays are removed using LACosmic (van 
Dokkum 2001). For each arm, we extract the orders and rec¬ 
tify them in wavelength space using a wavelength solution 
previously obtained from the calibration frames. The result¬ 
ing rectified orders are then shifted and co-added to obtain 
the final two-dimensional (2D) spectrum. We extract a one¬ 
dimensional spectrum (ID) from the resulting 2D merged 
spectrum, using our own IDL code that also produces the 
corresponding error file and bad pixel map. The final signal- 
to-noise ratio in the LTVB+VIS spectrum varies from sys¬ 
tem to system, with a minimum of ~ 50 and a maximum of 
~ 145 per Angstrom (0.2 A/pix). We finally perform an it¬ 
erative sigma-clipping to clean the spectrum of any residual 
bad pixels, sky lines and cosmic rays. No telluric correction 
has been applied so that prominent atmospheric absorption 
bands can still be seen in the final UVB+VIS spectra shown 
in rest-frame wavelengths in Figure 2(a) and Figure 2(b). 
For some of the systems, it is also possible to see emission 
lines from the background source, although the chosen PA 
minimizes their contribution to the stellar spectra. We mask 
these lines as well as telluric absorption when extracting the 
spatially-resolved kinematic profile. 

Because the near-infrared spectrum still suffers seri¬ 
ously from sky-line residuals that are difficult to remove us¬ 
ing the current pipeline, we limit ourselves in this paper to 
the UVB-VIS region of the spectrum and defer a full analysis 
of the infrared data to future publications. 

4 STELLAR KINEMATICS 
4.1 Global Kinematics 

We measure the luminosity-weighted line-of-sight velocity 
dispersion (LOSVD) of the lens galaxies from the final ID 
UVB-VIS spectra extracted from a rectangular aperture of 
2" x 1'.'5 centred on the galaxy. We use the Penalized Pixel 
Fitting (pPXF) code of Cappellari & Emsellem (2004) to 
determine the combination of stellar templates which, when 
convolved with an appropriate LOSVD, best reproduces the 
galaxy spectrum. The best-fitting parameters of the LOSVD 
are determined by minimizing a \ 2 penalty function, yield¬ 
ing the mean velocity and velocity dispersion (u and <r*, 
respectively) and their uncertainties. In Figure 3 we show 
examples of pPXF fits for the highest and lowest signal-to- 
noise systems. 

3 In the NIR arm, the total exposure time is the product of DIT 
(Detector Integration Time), NDIT (Number of DITs) and NINT 
(Number of INTegrations): T = DIT X NDIT X NINT. We used 
for each of the 3 scientific exposures 3 DITs of 278 sec. 
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Figure 2. (a). Final luminosity-weighted UVB-VIS ID X-shooter spectra (in black) and respective errors (in blue) extracted from a 
rectangular aperture of 2 " x 1//5 centered on the galaxies, shifted to the rest-frame and smoothed (3-pixel boxcar) for displaying 
purposes only. IDs, redshifts of the lens and average S/N (per A) are shown in the panels. Flux is in units of 10 _16 erg/s/cm 2 /A. Telluric 
absorption lines and emission coming from the background source have not been removed from the spectra. The peak in the error spectra 
between 4500 and 5500 A shows the point at which UVB and VIS ID spectra have been joined. 
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Figure 2. (b). continued 
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Table 1. Properties of observed systems. 


SLACS System 

Zlens 

zbg 

Reff(kpc) 

REin (kpc) 

m v (mag ) 4 

PA slit (°) 

SDSSJ0037-0942 

0.1955 

0.6322 

7.03 

4.95 

16.90 

11.4 

SDSSJ0044+0113 

0.1196 

0.1965 

5.56 

1.72 

16.32 

151.3 

SDSSJ0216-0813 

0.3317 

0.5235 

12.6 

5.53 

18.36 

81.2 

SDSSJ0912+0029 

0.1642 

0.3239 

10.8 

4.58 

16.56 

11.7 

SDSSJ0935-0003 

0.3475 

0.4670 

20.7 

4.26 

17.71 

145.2 

SDSSJ0936+0913 

0.1897 

0.5880 

6.61 

3.45 

17.12 

145.3 

SDSSJ0946+1006 

0.2219 

0.6085 

8.33 

4.95 

17.78 

10.3 

SDSSJ1143-0144 

0.1060 

0.4019 

9.21 

3.27 

15.83 

118.7 

SDSSJ1627-0053 

0.2076 

0.5241 

6.66 

4.18 

16.91 

6.9 

SDSSJ2343-0030 

0.1810 

0.4630 

8.27 

4.62 

17.17 

144.0 

SDSSJ1148+1930 1 

0.4440 

2.3815 

12.5 

29.0 

20.02 

99.0 


1 This system is not selected from SLACS but it satisfies the LRG criteria and was observed in a 
similar set-up during the XLENS pilot program (Spiniello et al. 2011). 

2 From the F606W filter on WFPC2. 


We focus on the absorption lines between 3500-6500 A 
(including Ca K and H, G4300, H^, Mg;,, strong Fe lines, 
NaD and some TiO molecular absorption bands). To min¬ 
imize errors due to potential mismatches between the res¬ 
olution of the templates and the galaxy spectrum, we use 
X-shooter spectra obtained as part of the X-shooter Stel¬ 
lar Library (XSL) survey 4 (Chen et al. 2014), with higher 
instrumental resolution. For the galaxy spectrum we use a 
1.5-arcsecond slit (corresponding to (<7mstr) ~ 35 kms -1 , 
in the observed frame, and to (cri ns tr) ~ 28 kms -1 in the 
restframe) while for the stellar templates the slit width are 
0'.'5 and 0 , /7 in LTVB and VIS respectively, corresponding to 
(crinstr} ~ 12 kms -1 . As a test of the accuracy of our mea¬ 
surements, we use the more heterogeneous MILES 5 stellar 
template library by Sanchez-Blazquez et al. (2006). We se¬ 
lect 100 stars (F, G, K, M) in the range 3525-7500 A, with 
2.5 A FWHM spectral resolution (corresponding roughly to 
30 ^ ctmiles ^ 70 kms -1 in the selected wavelength range), 
which is comparable to the XLENS spectra resolution in 
the selected regions (ctxsh ~ 35 kms -1 ) that have been 
smoothed to the same spectral resolutions with a Gaussian. 
Fitting with different stellar templates does not appreciably 
affect the resulting LOSVD. As a second test, we fit the same 
stellar templates to two different spectral regions for each 
galaxy separately (blue region: 3800 — 5000 A and red region: 
5000 — 6500 A). Only for two systems (SDSSJ0936+0913 
and SDSSJ1143-0144) do we find slightly different results 
between the blue and the red parts of the spectrum, but 
the results are always consistent within 2cr. The scatter be¬ 
tween the different fits is used to estimate additional system¬ 
atic uncertainties related to template mismatch and spectral 
coverage. These errors are folded into the final error budget. 
We report our results and the comparison with other pub¬ 
lished velocity dispersion results for each ETG in Table 2. 
We find a good agreement between our values extracted from 
a rectangular aperture of 2" x l'/5 and the SDSS velocity- 
dispersion measurements obtained in the S^-diameter spec¬ 
troscopic fiber (for most of the systems within 1 <t error). We 
also find good agreement with the results by Czoske et al. 

4 http://xsl.u-strasbg.fr 

5 http://www.iac.es/proyecto/miles/pages/stellar- 
libraries/miles-library.php 


(2012) obtained from (a larger) aperture-integrated spectra 
from the VIMOS Integral Field data on five systems that 
overlap with the XLENS sample. 

4.2 Spatially-Resolved Kinematics 

To extract the spatially-resolved kinematic information, we 
define a number of spatially-varying apertures (with ade¬ 
quate S/N ratio, S/N > 10) along the radial direction for 
each galaxy, and we sum the signal within each aperture. 
Apertures are always defined to be larger than the seeing 
in order to have approximately independent kinematic mea¬ 
surements for each aperture. The streaming motion and ve¬ 
locity dispersion are measured from each spatially-resolved 
spectrum using pPXF as described above. Also in this case 
we perform tests to check our error determinations: we fit 
the blue and red spectral regions separately, masking out the 
most prominent telluric lines in the VIS-red range; and we 
use different stellar templates from the two different stellar 
libraries (eleven XSL stars of G, K and M spectral types and 
a sub-selection of MILES stars of the same spectral types). 
The uncertainties on the inferred kinematics are estimated 
by adding in quadrature the formal uncertainty given by 
pPXF and the scatter in the results for different templates 
and spectral regions. 

Our goal is to determine the kinematic profiles up to 
about one effective radius for all the systems. The first 
and second moments of the LOSVD out to their effective 
radii, complemented with two-integral axisymmetric models 
(Barnabe & Koopmans 2007; Barnabe et al. 2009) allow us 
to obtain a precise measurement of the logarithmic slope of 
the total mass density profile of the lens galaxy (e.g. Czoske 
et al. 2008; Barnabe et al. 2011, 2013; Spiniello et al. 2011). 

The rotation and the velocity dispersion profiles for the 
ten XLENS-SLACS galaxies are shown in Figure 4. The up¬ 
per line shows systems with moderate rotation, whereas the 
lower line shows systems with almost-flat velocity profiles. 
The weighted average values are always consistent within 
the formal error with the luminosity weighted values for an 
aperture of 2" x l'/5. 

We infer the averaged angular momentum (Ah) using 
the definition of Emsellem et al. (2011), with the spatially- 
resolved rotation velocity and velocity dispersion. We derive 
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Table 2. Luminosity-weighted stellar kinematics of the lenses, ( ctxsh ) has been extracted from a rectangular aperture of 2 " x l'/5 
centered on the galaxies, while ctsdss is measured within an aperture with a diameter of 3 " . VIMOS velocity dispersions were measured 
from aperture-integrated spectra from data from the VIMOS Integral-Field Unit (Czoske et al. 2012). 


ID name 

(^xsh) (kms -1 ) 

ctsdss (kms -1 ) 

0VIMOS (kms 1 ) 

SDSSJ0037-0942 

277 ± 6 

279 ± 14 

245.3+£;2 

SDSSJ0044+0113 

260 ±8 

266 ± 13 

SDSSJ0216-0813 

327 ± 19 

333 ± 23 

340.7+77 

SDSSJ0912-I-0029 

325 ± 10 

326 ± 16 

306.51 3 ?® 

SDSSJ0935-0003 

380 ± 22 

396 ± 35 

330.4+l^g 

SDSSJ0936+0913 

256 ± 18 

243 ± 12 

SDSSJ0946+1006 

300 ± 22 

263 ± 21 


SDSSJ1143-0144 

287 ± 18 

269 ± 10 


SDSSJ1627-0053 

303 ± 23 

290 ± 20 

272.6+g'g 

SDSSJ2343-0030 

298 ± 21 

269 ± 16 



4000 5000 6000 


Restframe Wavelength (A) 


J0216-0813 



4000 5000 6000 


Restframe Wavelength (A) 


Figure 3. Two examples of pPXF fits. The upper row shows one 
of the best cases (x 2 /DOF ~ 1.15), while the lower row shows 
the worst case (x 2 /DOF ~ 1.65). Top panels: UVB+VIS galaxy 
spectrum (black) and correspondent best-fit template (red), both 
in restframe wavelength. Bottom panels: Residuals from the fit. 
Bad-pixels excluded from the fitting procedure (sky line, telluric 
lines) are shown in blue. See the text for more information. 


the radial Xr profiles (left panel of Figure 5) and plot them 
as a function of ellipticity (e), calculated from the lensing 
models (e = 1 — q*, where q* is the axis ratio). This allows 
us to split the galaxy sample into slow- and fast-rotators. 
The A r parameter quantifies the (apparent) global dynam¬ 
ical state of a galaxy and has been used by the ATLAS 30 
team (Cappellari et al. 2011) to show that fast and slow ro¬ 
tators are physically distinct classes of galaxies with differ¬ 


ent orbital distribution and mass assembly histories. Simu¬ 
lations have demonstrated that galaxy mergers (major and 
minor) have a significant influence on the rotation prop¬ 
erties: major mergers produce a significant increase in the 
angular momentum of the dark matter haloes at large dis¬ 
tances from the center (Vitvitska et al. 2002) but it leads 
to a redistribution of the angular momentum of the central 
stellar component outwards (e.g. Bournaud, Combes & Jog 
2004). Fast rotators have either preserved or regained their 
specific angular momentum in the central part and there¬ 
fore they cannot have experienced major dry mergers which 
would expel most of the angular momentum outwards and 
avoid the building of a rotation-dominated disc-like struc¬ 
ture (Emsellem et al. 2007; Naab et al. 2014). Concluding, 
the A r, interpreted as a proxy of the specific angular mo¬ 
mentum, is an important quantity that allows one to link 
present-day kinematical profiles to the cosmological forma¬ 
tion history of a galaxy. 

The right panel of Figure 5 clearly shows that five of 
our galaxies are classified as slow rotators, whereas the other 
five are classified as fast rotators, but none of them show a 
very high A r value. It also appears that the most massive 
galaxies in our sample tend to rotate more slowly than the 
least massive ones (with the exception of J0912, which is 
classified as fast-rotator even though it has cr > 300 km s _ ). 
This result confirms the ATLAS 30 result that slow rotators 
dominate the high-mass end of ETCs. 


5 SUMMARY & CONCLUSIONS 

In this paper, we have presented our spectroscopic follow¬ 
up survey of massive early-type lens galaxies: the X-shooter 
Lens Survey (XLENS). In particular, we 

• highlight the scientific motivations and main goals of 
XLENS. 

• present the current sample of eleven massive early-type 
lens galaxies, discussed the selection criteria and described 
the data reduction process for the XLENS-SLACS galaxies. 
The remaining lens, also known as “The Cosmic Horseshoe”, 
has been used as a pilot program of the survey and is the 
subject of the first paper of the series (XLENS I, Spiniello 
et al. 2011. 

• present the spatially resolved kinematic profiles up to 
~ 0.8 — 1.5 .Reff for the lens galaxies. We found luminosity- 
weighted velocity dispersions mostly within ltr agreement of 
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Figure 4. Spatially-resolved kinematics of the lenses. Upper panels show the systems which have mild rotation and are classified as 
fast-rotators from Figure 5, lower panels show systems which do not rotate and are classified as slow-rotators. 


those calculated from SDSS. Our velocity dispersions also 
agree with the results of Czoske et al. (2012) extracted from 
VIMOS integral-field spectra on 5 systems in common. 

• derive the Xr spatial profiles and calculate \n e{{ /2 to 
split the sample into fast- and slow- rotators. We confirmed 
the ATLAS 3D results that slow rotators dominate the high- 
mass end of the galaxy mass function. 

The combination of high-S/N spectroscopy from X- 
shooter with a strong gravitational-lensing mass determi¬ 
nation enables us to conduct an in-depth study of the dark 
and luminous matter within ~ lR e ff and for the first time 
constrain the normalization, shape and cut-off mass (Mi ow ) 
of the low-mass end of the IMF. We make use of the state- 
of-the-art joint lensing and kinematic code (Barnabe et al. 
2012) and new SSP models (Conroy & van Dokkum 2012) 
that have been constructed specifically for the purpose of 
measuring the IMF slope down to ~ 0.1 Mq for old, metal- 
rich stellar populations. The version of the SSP models 
that we use have been extended beyond their original pa¬ 
rameter spaces, to explore super-solar [a/Fe] and super¬ 
solar metallicity and to disentangle elemental enhancements, 
metallicity changes and IMF variations in massive early-type 
galaxies (ETGs) with star formation histories different from 


the Solar neighborhood (following the method presented in 
Spinicllo, Trager & Koopmans 2015). 

In two following papers of the series, we use the joint 
lensing and dynamics code CAULDRON 6 to perform the 
analysis of the total mass distributions for the XLENS sys¬ 
tems. We put constraints on the internal dark matter frac¬ 
tions, the orbital structure, the shape and flattening of the 
dark matter halo and investigate possible correlations of 
these quantities with physical parameters of the lens sys¬ 
tems such as size, stellar density and stellar velocity dis¬ 
persion. Thanks to the high quality of our lensing and kine¬ 
matics data, this technique also makes it possible to put firm 
constraints on the stellar mass of the analyzed galaxies in¬ 
dependently of assumptions about their stellar initial mass 
functions (IMF) or knowledge of their stellar populations. 

Moreover, we apply to the sample the line-index-based 
stellar population analysis presented in Spiniello et al. (2014) 
to study in detail the stellar contents of the systems and 
to constrain their integrated ages, metal abundances and 

6 ‘Combined Algorithm for Unified Lensing and Dynamics Re- 
constructiON’, (Barnabe & Koopmans 2007; Barnabe et al. 2009, 
2012) 
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Figure 5. Left: Radial A/,■ profiles. The 1-sigma contours are showed with colored shaded regions. Right: The A/c rf /2 — e relation for 
the XLENS galaxies. The dotted black line shows the division between slow- (below the line) and fast (above the line) rotators within 
R a ff/ 2 - All most massive galaxies in our sample, except for J0912, belong to the slow-rotator group. 


IMF slopes and to search for possible correlations between 
stellar population parameters and structural properties of 
the systems (i.e., mass, size, stellar density and/or velocity 
dispersion). 

Finally, in a Letter already submitted, we study the 
impact of a non-universal IMF slope and/or a non-universal 
IMF low-mass cut-off on the inferred mass-to-light ratios 
by comparing the two completely independent constraints 
on the stellar mass, one from the fully self-consistent joint 
lensing+dynamics (L&D) analysis, and the other from spec¬ 
troscopy and SSP modelling. For this study we follow the 
approach presented in Barnabe et al. (2013) to obtain the 
joint 2D posterior PDF in the slope-Mi ow space. 

In the future, we plan to extend the current sample, ob¬ 
serving lenses within the same range of velocity dispersions 
(c*,sdss > 250 kms - ) at higher redshift (up to 2 ~ 1) 
to test possible evolution of the IMF with cosmic time. We 
also plan to observe low-mass ETGs to build a complemen¬ 
tary dataset to our more massive ETGs with the purpose 
of covering a wider stellar mass range and to test possible 
variation of the IMF slope with structural and stellar param¬ 
eters, other than the stellar mass of the galaxy (e.g. total 
dynamical density, as we recently suggested in Spiniello et 
al. 2015). 
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